Microscopic mechanisms of dephasing due to electron-electron interactions 
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We develop a non-perturbative numerical method to study tunneling of a single electron through 
an Aharonov-Bohm ring where several strongly interacting electrons are bound. Inelastic pro- 
cesses and spin- flip scattering are taken into account. The method is applied to study microscopic 
mechanisms of dephasing in a non-trivial model. We show that electron-electron interactions de- 
scribed by the Hubbard Hamiltonian lead to strong dephasing: the transmission probability at flux 
$ = n is high even at small interaction strength. In addition to inelastic scattering, we identify 
two energy conserving mechanisms of dephasing: symmetry-changing and spin-flip scattering. The 
many-electron state on the ring determines which of these mechanisms will be at play: transmitted 
current can occur either in elastic or inelastic channels, with or without changing the spin of the 
scattering electron. 

PACS numbers: 73.63.-b 72.10.-d 71.10.-w 
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I. INTRODUCTION 

Advances in the semiconductor technology made it 
possible to study quantum interference effects in meso- 
scopic systems where the wave nature of electrons plays 
an essential role. Particularly noteworthy are the studies 
of the Aharonov-Bohm (AB) oscillations in mesoscopic 
ringsii2i£. The analysis of results in terms of the single- 
electron picture turned out to be inadequate to describe 
the totality of phenomena. Inelastic scattering of elec- 
trons is believed to be the predominant mechanism re- 
sponsible for the loss of the phase coherence in such ex- 
periments and the suppression of the h/e oscillations. 
When an electron interacts with optical phonons, the 
dephasing only occurs through inelastic processes^. At 
low temperatures the phonon degrees of freedom freeze 
out, therefore other mechanisms for dephasing should 
be taken into account. Measurements of the dephasing 
time saturation at low temperatures^ show that zero- 
point fluctuations of the electromagnetic environment 7 
could play a role in explaining this anomalous behav- 
ior. It is nevertheless believed, that at low temperatures 
the electron-electron interaction is the dominant mech- 
anism for dephasing^. Further support for the impor- 
tance of electron correlations at very low temperatures 
comes from recent measurements of anomalous temper- 
ature dependence of the dephasing time in mesoscopic 
Kondo wiresifi where non- Fermi-liquid behavior has been 
found below the Kondo temperature. 

The AB geometries have been theoretically studied us- 
ing self-consistent mean-field approximations that break 
down for degenerate levels, which physically happens 
at very low temperaturesii*i2*ii The mean-field ap- 
proximation does not describe transitions in which the 
symmetry of the many-electron wavefunction of corre- 
lated bound electrons changes and it is thus inadequate 
to study dephasing. Kcldysh type Green's functions 
and numerical renormalisation group techniques'**^ or 
equation of motion methodic have been applied to AB 
systems, where calculations were limited to interacting 



quantum dot with two levels coupled to reservoirs. Par- 
ticular attention was devoted to the appearance of the 
Kondo physics induced by changing the magnetic flux, 
however no spin-flip induced dephasing has been investi- 
gated by these methods. A study of a Coulomb blockade 
regime was recently done by Xiong and Xiongi 7 - with a 
method, similar to the one proposed in the present work. 
Their Hamiltonian, however, maps on a non-interacting 
model in the limit when the coupling to the leads is 
zero. Furthermore, they have only investigated spin- 
less fermions and neglected inelastic scattering channels. 
Transmission of two interacting electrons was recently 
studied on the basis of continuous two-particle Hamilto- 
nian where an enhancement of transmission with increas- 
ing interaction strength was founcU&. 

To shed some new light on the problem of dephasing 
in electron tunneling experiments, there is obviously a 
demand for a capable method that would treat the prob- 
lem of the scattering of an electron through a finite re- 
gion where electron-electron (e-e) interactions would be 
exactly taken into account. Such a method should be 
based on the use of exact correlated many-electron wave- 
functions. 

In this Paper we propose a method that treats the e-e 
interactions by direct diagonalisation of the many-body 
Hamiltonian using iterative (Lanczos) technique. The 
method naturally takes into account spin-flip processes, 
so it can predict the ratio of spin-flip over normal scat- 
tering processes. This makes the technique interesting 
for calculating spin-polarized transport^ in the field of 
spintronics. 

We apply the method to study single-electron trans- 
mission through an Aharonov-Bohm (AB) ring with e-e 
interactions. As it is widely known, an electron perfectly 
reflects from an AB ring when the flux $ penetrating the 
ring is such that the phases gained by the electron trav- 
eling through the lower or the upper arm of the AB ring 
cancel each other (<fr = 7r). Such reflection occurs for any 
energy of the incident electron. This remains true even 
when there are electrons bound on the AB ring, as long as 



the system remains non-interacting. The main purpose of 
our investigation is in the influence of the finite Coulomb 
repulsion on the transmission of the electron in the case 
described above. We choose the Hubbard model to de- 
scribe the AB ring. The Hubbard model is the simplest 
yet most important nontrivial prototype model for corre- 
lated electrons in the solid state. As we will show, finite 
Coulomb interaction in certain cases leads to finite trans- 
mission of the incident electron, despite the fact that the 
total wavefunction for the scattering electron and elec- 
trons bound on the AB ring preserves full quantum coher- 
ence. We will therefore refer to the processes that lead to 
finite transmission in the case when <E> = tt as dephasing 
processes since they clearly lead to diminished AB oscil- 
lations observed in experiments. We show that dephasing 
can occur either by a) inelastic processes where the tun- 
neling electron excites bound electrons on the ring or by 
b) elastic (in regard to energy) processes, where the tun- 
neling electron changes the symmetry or the spin of the 
degenerate many-electron wavefunction. No exchange of 
energy is required in the latter cas o 1 ^ 2 ?' 21 ! 22 : dephasing 
occurs because the tunneling electron leaves a trace on 
its "environment" , which consists of bound electrons. 

II. METHOD 

The proposed method is based on the multichannel 
scattering technique that was developed for studying the 
tunneling of a single electron in the presence of scattering 
by phonon a 23 i 24 . Since its introduction, it has been suc- 
cessfully applied to a variety of problems where a single 
electron is coupled to phonon mode a 4 ' 2 ^ 2 9' 2 7i 2 9i 2 ?i? l i ) i? 1 . 
It was even incorporated into Landauer theory where 
the influence of electron-phonon scattering on the non- 
equilibrium electron distribution has been investigated 32 . 
We now generalise this method to study many-electron 
problems. 

A. Model Hamiltonian 

While the method can be applied to more general sit- 
uations and arbitrary geometries of the interaction re- 
gion, we choose for simplicity a particular physical sys- 
tem which will also serve as a toy-model for the study 
of the e-e interaction induced dephasing. We thus con- 
sider an AB ring coupled to two ideal one-dimensional 
leads, see Fig. ^ The ring is described by a Hubbard- 
type Hamiltonian 

Hr'mg ^ Cj,a^ te ^ Cj-\-l,a^ Cj,t7 "T" h.C.) 

(1) 
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The operator creates an electron with spin a at 

site j and we make a formal identification = Cx >a K 
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Extended interaction region 




Left lead Aharonov-Bohm ring Right lead 



FIG. 1: Aharonov-Bohm ring. Magnetic flux penetrates the 
center of the ring. 

The phases <frj describe phase changes due to magnetic 
flux penetrating the ring. We chose a gauge in which we 
ascribe the total phase change due to magnetic field flux, 
$ = 27t<I>m/ < I > o, (where $m is the magnet field flux and 
$0 = h/e is the flux quantum) to a single element, e. g. 
4>i = $ and <j>j = 0, j ^ 1. 

The leads are described by a tight-binding Hamiltonian 

H\ ca d = -head ^ ai+l,ff^O»,(7 + n - c - (2) 

-tlead X! *b it<7 + h.C. (3) 

i,a 

The operator creates an electron with spin a at site 
i on the left lead, while the operator bi „ does the same 
on the right lead. The ring is coupled to the electrodes 
with coupling constants to: 

H c = -to (ai.o-^ci.CT + h.c.) - 1 ^2 {bi^Ci.a + h.c.J 

(X (T 

(4) 

The coupling constant to need not be small: the method 
applies equally well for strong coupling between the in- 
teraction region and the leads. 

B. Correlated many-electron states and exact 
diagonalisation 

The transmittivity can be meaningfully defined in 
many-electron scattering problem only if one single elec- 
tron leaves the scattering region. For this reason, we re- 
strict the energy of the incoming electron to be below the 
ionisation threshold. Our approximation then consists of 
taking into account only those many-electron states in 
which at most one (scattering) electron is located out- 
side the ring. Before the impact of the electron (which, 
for convenience, will be chosen to have spin up), there 
are n = ri| + n± other electrons bound on the AB ring. 
We truncate all many-body states, where additional elec- 
trons hop from the interacting region to the lead. When 
physical parameters of the system, e.g. (e, t, U), are cho- 
sen in such a way that these n electrons are bound in the 
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interacting region, the approximation is equivalent to ne- 
glecting the exponentially decaying tails of the n-clcctron 
wavefunction in the leads. 

Before the scattering, the bound electrons are therefore 
in one of the n-particle eigenstates of the Hamiltonian 
-ffring, Eq. QJ. We denote these states by \a\) and their 
energies by ej. The superscript index | denotes that the 
electron in the lead has spin up. 

When the incoming electron enters the ring, the sys- 
tem is in a superposition of the n+ 1-particle states which 
we denote by \(3i). These states are not necessarily eigen- 
states of -ffring- After the scattering there is a single elec- 
tron in the leads, while the ring is in a superposition of 
the n-particle eigenstates of H rmg . These states are the 
\cq) states and (if the spin of the scattering electron can 
be flipped) the n-electron eigenstates with rij + 1 spin-up 
electrons and nj_ — 1 spin-down electrons. These spin- 
flipped states are labeled by |q|) and their energies by ej. 
Because all the possible states of the ring after scattering 
are orthogonal to each other, the outgoing channels are 
well defined and the current is conserved. 

We calculate the eigenstates \af) using exact diago- 
nalisation of the Hamitonian H vlVLg in the suitable region 
of the many-particle Hilbert space, taking into account 
that the Hubbard Hamiltonian is invariant with respect 
to rotations in the spin space. The diagonalisation is 
therefore performed in the constant (n, S z ) space, where 
S z is the conserved component of the total spin in the 
z direction. The method can be applied to Hamiltoni- 
ans that don't conserve S z at the expense of significantly 
more time consuming numerical calculations. 

At zero temperature, the electron scatters on the 
ground state of the n-particle state in the ring, |aj). 
During the scattering the electron can loose energy by 
exciting the bound electrons into one of the excited \af) 
states. The probability of such transitions is a rapidly 
decreasing function of the energy loss, therefore only a 
small number of the scattering channels (states \af}) has 
to be considered. This observation is essential for numer- 
ical performance of the method: we can efficiently calcu- 
late the states from the bottom of the spectrum of the 
matrix representations of if in suitable (n, S z ) subspaces 
using the iterative Lanczos technique. We have used the 
implicitly restarted Lanczos method, as implemented in 
ARPACK package^. The eigenvalues and eigenvectors 
were computed to machine precision. 

By taking into account only the allowed states, the 
wave-function that describes the scattering of one elec- 
tron on the AB ring is given by 

oo oo 

i=l j,a i=l j,a /g\ 



where dfj a , dfj a and e k are the coefficients to be deter- 
mined. 



C. Reduction to a sparse system of linear equations 

We consider a steady-state scattering described by the 
the Schrodinger equation 



UW) = EW) 



(6) 



with H = Bring -\~H; i ea d -\-H c . Strictly speaking, this equa- 
tion cannot be solved in the restricted space, spanned on 
the \af) and |/3j) states, because applying the Hamil- 
tonian to the wave-function ansatz takes us out of this 
space by generating terms where more than one electron 
can be found outside the scattering region. Omission of 
these terms represents the main approximation used in 
our method. This approximation leads to an error which 
is not significant for suitably chosen model parameters 
(see below). 

Operating on the equation from the left with (/3;| 
we obtain 

where &'s denote scalar products 



(8) 



while hi t k — {Pi \H V i ng \Pk) are the matrix elements of 
Hamiltonian H r i ns in the n + 1 electron subspace. 
By operating with (aJ|ai j(T from the left we get 

-*w4,i, CT - to E( 6 L^)* e * + £ i rf L> = Ed lj,*- ( 9 ) 



D. Pruning the leads 



In an open outgoing channel (j, a) a plane wave can 
propagate, so that d\ ■ a = exp(ikj l0 -)d\ j a . By energy 
conservation the wave number kj jlT is obtained from 



e - 2ii cad cos(K) = ej - 2<i cad cos(fc 3i(T ). 



(10) 



The energy eo is the initial energy of the n-electron bound 
state on the ring, K is the wave number of the incoming 
electron, and is the final energy of the bound electrons. 
Equation @ can thus be written as 



E - e? + ticad exp(ikj 



(11) 



Similar equation can be obtained for exponentially de- 
caying (closed) outgoing channels that we also take into 
account (up to some cut-off energy, above which the in- 
clusion of further closed channels does not change the 
results). These are defined through the relation d\ ■ a = 



cxp(-K J> )c?5 J j a and 

e - 2tiead cos(iT) =e°- - 2ii cad cosh(kj 



(12) 
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Equation © can in this case be written as 



E - 



+ t\ ea d exp( 



(13) 



In the incoming channel we have both the incoming and 
outgoing waves, Q * = exp(— iKm) + r exp(iKm). We 

obtain cS^ot = exp(ii^)ci^ ^+exp(— UK) — 1. The equa- 
tion for the incoming channel thus contains an additional 
inhomogeneous term exp(— UK) — 1 and Equation for 
the incoming channel is 



r1 L 

i,o, T 



_ *o Y,k( b k,o,i)* e k ~ *ioad[exp(-2ii4:) - 1] 
E-e Q 



tiead exp(iisT) 

(14) 

Using Eqs. (|TT]) . l|13[l. ((141) . an d similar equations for 
the right lead, both leads can be removed (pruned) from 
the problem^. 

Equations J7J|, l|ll|) and equivalent equations for other 
outgoing channels form a system of linear equations for 
unknowns d[ J - (T , • CT and e^. This sparse system is 
solved for different energies of the incoming electron using 
the SuperLU library^. 

The partial transmittivity through channel (J, a) is 
given by 



TjAE) 



sin(fej, CT ) , R |2 



sin(2<0 



(15) 



Since the method is based on exact solution of many- 
electron problem, we can compute transmission at arbi- 
trarily large values of U. 



E. Extended interaction region 

Results can be improved by extending the interaction 
region which is solved exactly by the Lanczos method by 
adding additional sites from the leads. This procedure 
takes into account the decaying tails of bound electron 
wavefunctions in the leads at the expense of increasing 
the computational Hilbert space. The error due to the 
omission of the terms corresponding to a second electron 
jumping out of the original interaction region [see discus- 
sion following Eq. ©] is exponentially reduced with the 
inclusion of each additional site from the leads. 

These improvements mainly lead to energy shifts of the 
resonance peaks while the general characteristics of the 
spectra remain unchanged. In principle, the region could 
be extended until the desired convergence is achieved. In 
our calculations the interacting region consisted of the 
ring and one additional site from each lead, see Fig. ^ 
In fact, we had to take into account the site on the left 
electrode in order to ensure that the incident electron 
spreads into two partial waves that propagate through 
both arms of the ring. The additional site on the right 
electrode is required so that the partial waves can in- 
terfere, which leads to the Aharonov-Bohm effect. The 



inclusion of these two sites was therefore essential in our 
studies of the dephasing mechanisms. 

In cases where the ground state of the interaction re- 
gion was degenerate, we averaged the transmittivity spec- 
tra over all the degenerate states. The variational space 
taken into account in our calculation was equivalent to a 
Hubbard model on 8 sites with no translational symme- 
try. 

The largest problem that could be effectively solved 
has 3 bound electrons with spin up and 4 electrons with 
spin down. In this case there are ~8000 states \af) and 
-5000 \(3i) states. We kept 200 lowest lying \af) states 
to define our scattering channels (diagonalization took 
3 minutes of a modest personal computer). Solving the 
sparse system of complex linear equations for a range of 
the incoming electron energies took 200 minutes for 237 
data points (or about a minute per data point on the av- 
erage) . This step is the most computationally demanding 
part of our technique. This is the main reason why we 
limited our calculations to 8 site models, even though the 
Lanczos method easily handles much larger lattices. 



III. RESULTS - ONE BOUND ELECTRON 

We now investigate the effect of interactions on an 
electron as it tunnels through the ring. The incoming 
electron has spin up and there is one bound electron 
with spin down inside the ring. The on-site energies are 
e = — 4.5iicad, the overlap integrals are t = v^icad, and 
we set iicad = 1. 

First we consider the non-interacting case. In the ab- 
sence of the magnetic field the transmission reaches unity 
at the resonance, Fig. The electron is fully reflected 
at any incident energy when the magnetic flux is <I> = n, 
Fig. 03. This is the usual Aharonov-Bohm effect. 
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FIG. 2: Transmission probability as a function of the incident 
electron energy for one electron with spin down bound on the 
ring. The incoming electron had spin up. The coupling to the 
lead is to = 0.4. In all cases transmission is purely elastic. 



5 



We now turn on the interaction. At <f> = we still 
see a unity peak at the energy of the single-electron res- 
onance, followed by smaller satellite peaks caused by the 
interaction, Fig. At <f> = ir, when in the absence 
of U the electron is fully reflected, we obtain very high 
transmission probability despite relatively small U = 0.2, 
Fig. 01. In the largest peak the transmission approaches 
the value T = 1/2. Since the incoming electron and the 
bound electron are not entangled, their total spin is not 
well defined, therefore the total wavefunction is a super- 
position of a singlet and a triplet state with equal ampli- 
tudes: |U) = 1/V2(|5 = 1,5, =0) + \S = 0, S z = 0)). 
The triplet scattering has zero transmission probability 
at $ = 7r since in the Hubbard model only singlet elec- 
trons interact. The singlet scattering, however, reaches 
the unitary limit at the main resonance peak. Averaging 
over both contributions, we indeed get T = 1/2. 

The spin-flip scattering part of the transmission prob- 
ability is shown in the inset in Fig. Eli- The spin-flip and 
normal scattering contribute equally to the total trans- 
mission probability. Both are purely elastic with respect 
to energy changes. 



form two linear combinations of states: 
1 




e nJ + e J e m c 



J,<r> 



(18) 




3,v 



It is easy to see that the coefficient of c m / 2 , CT m the ex- 
pression for aL^ nt<T is zero and likewise for the coefficient 
of co,(T in the expression for a,R. nM . This means that 
the eigenstate denoted by L is coupled only to the left 
electrode, while eigenstate R is coupled only to the right 
electrode, see Fig. [21 In the non-interacting case the in- 
coming electron can only tunnel from the left electrode 
to an L state. This state is decoupled from the right elec- 
trode and since there is no term in the Hamiltonian, that 
would allow transitions from L to R state, the electron 
is fully reflected. 



Lead decoupling at $ = 7r and scattering 
mechanisms at non-zero U 



To gain more insight in the mechanism of non-zero elec- 
tron tunneling probability, we present a simple physical 
picture of electron tunneling for the case of $ = it. We 
first transform the A-B ring Hamiltonian from the 
basis of localized states to a basis in k-space. For more 
generality, we can now assume that the ring consists of 
an arbitrary even number to of sites, which we now num- 
ber from to to — 1, so that the site is coupled to the 
left electrode, while the site to/2 is coupled to the right 
electrode. 

The non-interacting part of the Hamiltonian is di- 
agonal in the plane-wave basis, 



(16) 



3=1 



with wave numbers k n given by the periodic boundary 
condition exp(i/c„TO,) = 1, or k n — 2im/m, where n = 
0, ±1, . . . , ±to/2— 1, m/2. The corresponding eigenvalues 
are 



E n = e — 2icos(fc„ — $/to) 



(17) 



When $ = 7r, all the non-interacting eigenstates are 
twofold degenerate since cos(fc„ — Z-) = cos(/ci_ n — 
The complete orthonormal set of states is therefore com- 
posed of to/2 pairs of states with wave-numbers k n and 
fci_„ for n ranging from 1 to to/2. For each pair we can 




FIG. 3: At $ = 7T, one can rotate each pair of the degenerate 
eigenstates in such a way that one of them couples only to the 
left electrode and the other to the right electrode. In absence 
of interactions both electrodes are effectively decoupled. 

There are therefore two equivalent physical descrip- 
tions of zero transmittivity of an AB ring. One can ei- 
ther consider it as a destructive interference of partial 
electron waves that travel in the upper and lower arm of 
the ring, or as an effective decoupling of both electrodes 
due to a topological phase shift. 

We will now write the interacting part of the Hamil- 
tonian in the new basis and search for processes that are 
responsible for the nonzero transmission. From Eqs. 
and l |18| l we can express Cj >tr as 



2e i(^) m / 2 , 

c,- jCr = — V" ( cos(fc„j) a i.«.cr + s\n{k n j)a R . ni<J 
\>2m — ' v 

v n— 1 



(19) 
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where we have introduced shifted wave numbers 

kn = -(n-l)- (20) 
The particle number operator can then be expressed 

as 

2 

c j,<T C 3,a = — ^ cos C k P 3) c os(k q j)a f Lp!7 a L ^ a 
p,q 

+ sin(fcpj) sm(k q j)a Rp ^a R , q , a ( 2 1) 
+ cos(k p j) sax(k q j)a\ p a Rtqta 
+ sm(k p j) cos(k q j)a Rp an, g ,<T- 

We now see that the Hubbard interaction term 

Y^j c j,V c j,1 c j,l c j,l ^ s a sum over 3 01 products of four 
trigonometric functions. Each one of these products can 
be written as a sum of trigonometric functions by using 
trigonometric reduction formulas such as, for example, 

8 sin(a) sin(6) sin(c) sin(d) = 

— cos(a — b — c — d) + cos(a + b — c — d) 

+ cos(a — b + c — d) — cos(a + b + c — d) (22) 
+ cos(a — b — c + d) — cos(a + b — c + d) 

— cos(a — b + c + d) + cos(a + b + c + d) . 

We note that reduction formulas for an even number of 
sine and cosine functions consist of a sum of cosine func- 
tions, while the reduction formulas for an odd number of 
sine and cosine functions consist of a sum of sine func- 
tions. This fact is important to understand the selection 
rules that lead to dephasing. 

The arguments of functions on the right hand 
sides of the reduction formulas are sums of the form 

^(p-!±(9-!)±(r-!)±(s-!)), Le - of form 

-22^ ( where t is an integer. When the summation over 
site index j is performed most of the terms will drop, 
since 




All interaction terms with a coefficient that after 
trigonometric reduction involves a sine function will 
therefore vanish. Such vanishing terms come from prod- 
ucts of an odd number of trigonometric functions of each 
kind, therefore they are of the form such as 

a R,pA a L,q,ta L ,r,l a L,8,l- (24) 

Such terms would allow (for p — q, r — s) transitions of 
the tunneling electron from state L to R without chang- 
ing the bound electron state (i.e. without leaving any 
imprint on the environment). Such transitions would 



clearly be in contradiction with our understanding of the 
dephasing in the A-B rings. 

The terms with four sine and with four cosine func- 
tions are of little importance for our purposes. They de- 
scribe intcrlevcl repulsion and intcrlcvcl transitions with- 
out changes of the L/R character of the electron states 
and therefore do not lead to a finite transmission. We 
will focus instead on terms with two sine and two cosine 
functions. They are of three kinds. The first one consists 
of terms of the form 

a R,p,1 a L,g,ia{,r : i a R.sd- (25) 

These terms describe what we call symmetry-changing 
transitions: the tunneling electron (with spin up) in the 
L state jumps to a R state, while a bound electron un- 
dergoes a transition from RtoL state. Such a transition 
can either be elastic (with respect to the energy of the 
tunneling electron) if p = q, r = s (Fig. UJ, or inelastic 
(Fig. 13 . 

,---„L R L R 




a) before b) after 



FIG. 4: Transmission due to elastic symmetry-changing scat- 
tering. 



,---^L R L R 




a) before b) after 



FIG. 5: Transmission due to inelastic symmetry-changing 
scattering. 

The second kind of terms is of the form 

a L P ,T a £^T a kr,i a «. s <i- ( 26 ) 

These terms correspond to spin-Hip transitions: the tun- 
neling electron with spin-up in the L, q state makes a 
transition to a lower laying state L, p state, while a bound 
electron undergoes a transition from the R, s to the i?, r 
state, Fig. |SJ Transition can again be either clastic (with 
respect to energy) for p = q, r = s (Fig. |SJ, or inelastic 

(Fig- 13- 

Finally, terms of the form 

a R, P A a L,qAa R r A a L ,sA (27) 
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a) before b) after 

FIG. 6: Transmission due to elastic spin-flip scattering. 

, — . L R L R 



a) before 



b) after 



a) before b) after 

FIG. 7: Transmission due to inelastic spin-flip scattering. 

can correspond either to symmetry-changing (Fig. |SJ or 
to spin-flip transitions (Fig.[5J), depending on the p, q, r, s 
quantum numbers. 

L R L R 



a) before 



b) after 



FIG. 8: Transmission due to elastic symmetry-changing scat- 
tering (of the second kind). 



B. Scattering of a wave packet 

To illustrate more in detail our results presented in 
the previous subsection, we present here numerically ex- 
act calculation of scattering of an incoming electron (de- 
scribed as a wave packet with a given finite energy width) 
on the electron with the opposite spin, bound on the 
Aharonov-Bohm ring. Since we are dealing with a sim- 
ple case of only two electrons, this problem can be solved 
numerically exactly by direct integration of the two-body 
Schrodinger equation ihd\i/j) / dt = H\ip). We take into 
consideration a sufficiently high number of chain sites, so 
that the positional spread of the wave packet is smaller 
than the length of the left and the right lead. We chose 
N = 200 sites, where the 6-site Aharonov-Bohm ring oc- 
cupies positions ranging from 101 to 106. 

We construct the wavefunction at the initial time as 



4, 



JO). The operator 



fct 



creates 



y TpackotVTboundl"/- ^ Vj.bound 

an electron with spin down in the bound eigenstate of 



FIG. 9: Transmission due to elastic spin-flip scattering (of the 
second kind). 



the Aharonov-Bohm ring. We calculated this state using 



direct diagonalisation. The operator ip 



t 

f packet 



t 

1 packet 



C^exp 



(k - fc ) 2 



2a 2 



exp(- 



(28) 

where c| fe = 1/VN Ej=i Gxp(ikj)cy arl d C is a normal- 
ization constant. This operator creates an electron with 
spin up in a wave packet centered at site iV cen t e r which 
has the average wave number ko and a spread of a in the 
k-space. We choose fco = tt/2 to place the wave packet in 
the middle of the energy band of the leads with the group 
velocity v = dE/dk(k = k ) = t = 1. We set a = 0.13 

and iVcenter = 50. 

The equation of motion was then integrated using 
Bulirsch-Stoer method, which gives highly accurate re- 
sults for this type of problems. The accuracy and stabil- 
ity can be conveniently estimated by monitoring the devi- 
ation from the proper normalization of the wave-function. 
Using the Bulirsch-Stoer method the normalization dif- 
fers from 1 at the eighth decimal place after the scatter- 
ing. 

We set the parameters to e = —3.0, t = v3, to = 0.6, 
^lcad — !• For the non-interacting system the transmit- 
tivity at <& = it is for all electron energies, while the 
transmittivity of an interacting system with U = 1 is 
shown in Fig. The location and the spread of the 
energies of the wavepacket are represented in the figure 
by a two-sided arrow. 

The electron density before and after the scattering at 
<f> = 7r is shown in Fig. El for the non-interacting case 
and in Fig. ^| for the interacting case^t For [7 = 0, 
the wave-packet is perfectly reflected, as expected. For 
U = 1 the wave-packet is partially transmitted through 
the scattering region. In fact, the expectation value to 
find an electron in the second electrode, Pr, corresponds 
to the following average: 



Pe 



dkT(e(k))\^ packct (k,t = 0)\ 



(29) 



where the transmission T is calculated using the method 
from section II and is presented in Fig.^| This equation 
connects and thus validates the two distinct methods. It 
is furthermore worth stressing, that the probability of 
finding electrons with either orientation of spin in the sec- 
ond electrode is equal, see Fig. ^21 This can be explained 
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FIG. 10: Transmission probability as a function of the inci- FIG. 12: Electron density before and after the scattering of 
dent electron energy for one electron with spin down bound the wave-packet on Aharonov-Bohm ring at $ = ir: interact- 
on the ring. ing case. 



as follows: finite transmission is a direct consequence of 
interaction U which in the case of two electrons acts only 
on the singlet part of the wavefunction. The triplet part 
does not feel U due to the on-site nature of the inter- 
action. Transmission therefore occurs only through the 
singlet channel. 



II. OS 
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FIG. 11: Electron density before and after the scattering of 
the wave-packet on an Aharonov-Bohm ring at <E> = ir: non- 
interacting case. Note that the vertical scale is the same for 
both spin projections: the scale was chosen so that the wave- 
packet is clearly visible. 



C. Aharonov-Bohm oscillations 

Aharonov-Bohm effect is experimentally observed as 
magnetic flux dependent oscillations of the electric cur- 
rent through a mesoscopic ring structure*. From cal- 
culated T(E) spectra we could estimate the zero-bias 
conductance as G — GqT(Ep), where Go = 2e 2 /h is 
the conductance quantum and Ep is the common Fermi 
level of both leads. In our minimal model with a dis- 
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FIG. 13: Time dependence (for interacting case) of the prob- 
ability to find an electron with given spin projection either in 
the ring, in the left electrode or in the right electrode. 



Crete number of resonance states in the ring, the energy 
shifts of the peaks when the flux is changed (see Eq. iJHJl) 
lead to pronounced conductance variations not necessar- 
ily connected to the Aharonov-Bohm effect itself. This 
is a direct consequence of using a small discrete num- 
ber of sites. It is therefore more revealing to observe the 
variations of the integral of transmittivity over the whole 
energy band, J T(E)dE. This quantity is relatively in- 
sensitive to energy shifts of the peaks, while it should 
clearly show A-B oscillations which affect the height of 
all of the peaks. 

In Fig. [21 we present this integral as a function of 
$ for a number of increasing interaction strengths U — 
0, 0.1, . . . , 1.0. The amplitude of A-B oscillations notice- 
ably decreases as the interaction grows stronger. Fig- 
ure ^] also shows, that the integral transmittivity is es- 
sentially interaction independent around zero flux, $ = 0. 
A similar insensitivity of the transmittivity sum-rule has 
been discovered in the case of tunneling in the presence of 
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FIG. 14: The flux dependence of integrated transmission 
probability for different interaction strengths. The coupling 
to the leads is to = 0.6, all other parameters are as before. 
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FIG. 15: Transmission probability as a function of the in- 
cident electron energy for raj (nj) electrons with spin up 
(down). Parameters are e = —4.5, t = V3, U = 1.0, to = 0.3 
and t lcad = 1. 

electron-phonon coupling^Si3&. This insensitivity breaks 
down at larger U. 



IV. MANY BOUND ELECTRONS ON THE 
RING 

We now consider several interacting (U = 1) bound 
electrons on the ring. All presented cases are calculated 
at the flux value = 7T, unless otherwise specified. Spin 
of the incoming electron is up. We have limited the en- 
ergy of the incident electron to a half of the bandwidth, 
i.e. E — [—2,0], in order to avoid ionisation. Our main 
goal in this section is to investigate the circumstances, 
under which a scattering electron obtains a finite trans- 



mission probability at $ = n when scattering through 
the AB ring in the presence of many bound electrons. We 
show that in most cases Coulomb interaction leads to fi- 
nite transmission. In our work we refer to processes which 
cause finite transmission as dephasing processes. To avoid 
confusion we point out once more that the total wave- 
function describing a many-body state of the scattering 
electron and bound electrons preserves its full quantum 
coherence throughout the calculation. Our Hamiltonian 
does not contain coupling to external degrees of freedom 
that would naturally lead to dephasing. 

When the bound state on the ring consists of three 
electrons with spin up and one electron with spin down 
(Fig. 115b ) no spin-flip scattering is possible because such 
processes turn out to be energetically impossible. The 
ground state is however fourfold degenerate and the tun- 
neling electron can get through the ring at finite U by 
changing the symmetry of the many-electron state on the 
ring. Since the ground state is degenerate, this process 
is purely elastic. 

In the case of nj = 0, nj_ —4 (Fig. \15b ) the ground 
state is non-degenerate, however the spin-flip processes 
are energetically allowed. We therefore obtain transmis- 
sion probability only in spin-flipped channels. Since in 
this case the ground state is not degenerate, the trans- 
mission consists of purely inelastic processes. 

In the case when the ground state is degenerate and 
the spin-flip processes are allowed, we expect dephasing 
to occur both with or without spin flip. Such is the case 
of nj = l,ni = 3, Fig. 115b . The transmittivity without 
spin-flip is purely elastic, while the spin-flip processes 
are predominantly elastic, with small contribution from 
inelastic channels. 

Finally, for nj = 2,?ij — 2 electrons are fully reflected 
from the ring since there are no allowed scattering chan- 
nels in the appropriate energy interval, Fig. I15H . 

We finally show the influence of large U — 15 on the 
case of nj = 2, = 2, where at U = 1 transmission re- 
mained zero in the whole interval of incoming electron en- 
ergy due to widely spaced many-electron levels. At large 
[7 = 15 the energy difference between the non-degenerate 
ground state and the first excited state decreases in com- 
parison with U = 1 the states become com- 
pressed in the lower Hubbard band. We changed the 
on-site energy to e — —20 in order to keep the electrons 
bound on the ring. At $ = there are several energies 
at which the electron can resonantly tunnel through the 
ring, Fig. 116b . At $ = 7r, the electron can only tunnel 
inelastically. The energy difference to the first excited 
state in the n electron Hubbard band is approximately 
1.4. We find indeed that only the electrons that are more 
than 1.4 above the bottom of the energy band can tun- 
nel, Fig. 116b. Such inelastic processes occur both without 
(Fig. I16bl or with spin-flip (Fig. USUI. 
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FIG. 16: Transmission probability as a function of incident 
electron energy for n T = 2, = 2, (7 = 15 and e = —20. 



V. CONCLUSIONS 

Using a simple model and a new numerical method 
we have investigated physics of single electron tunnel- 
ing through the AB ring in the presence of correlated 
bound electrons. In particular we have focused on the 
the role of electron-electron interactions on dephasing. 
While the proposed method clearly has some limitations 
(small interacting regions, inability to describe ionization 
processes, neglect of many-body effects in the leads), it 
nevertheless allows to treat the strong-interaction prob- 
lem exactly and to identify the two principal microscopic 
mechanisms which lead to the loss of phase coherence 
in quantum interference experiments. We showed that 
a particle can tunnel through AB ring at $ = it elasti- 
cally by a) changing the symmetry of the many-electron 
state which is possible in the case of degeneracy or b) by 
flipping the spin. Tunneling can also occur in the inelas- 
tic channel by exciting the many-electron state on the 
ring into an excited state with or without the spin-flip. 



Depending on the number of bound electrons, their total 
spin, degeneracy of the ground state and available energy 
of the incoming electron, the total transmission can be 
composed of partial transmissions caused by either one 
of the listed processes. 

Using the novel method we have thus unraveled micro- 
scopic mechanisms based on electron-electron interaction 
that in a mesoscopic system contribute to a finite trans- 
mission through the AB ring in the case of $ = it. How- 
ever, since our method is based on small physical systems 
that can handle only a few lattice sites and interacting 
electrons we have no means at this stage to perform ac- 
curate calculation of the dephasing rate. 

Even though all presented results are obtained on the 
basis of zero-temperature calculations, the method can 
be generalised to finite temperatures with some addi- 
tional numerical effort. On the other hand, our results 
do not necessarily predict a finite dephasing rate at zero 
temperature. Since we treat only a single electron in the 
leads we are completely neglecting the effects of many- 
body interactions spreading from the interacting region 
to the electrons in the leads. This spread forms the basis 
for the Kondo effect. At temperatures below the Kondo 
temperature Tk, our approach therefore breaks down; in 
the Kondo regime the spins of the electrons from the in- 
teracting region couple into singlets with the electrons 
from the leads. This process prevents spin-flip scatter- 
ing, which in our calculation represents one of the mech- 
anisms for dephasing. Kondo coupling may also lift the 
degeneracy of the many-electron states in the interacting 
region and thus prevent transmission through the elastic 
channel which leads to dephasing at zero temperature 
according to our findings. Other mechanisms leading 
to dephasing in our approach might as well be modified 
in this low-temperature regime. We therefore conclude 
that despite the zero-temperature formalism used in our 
method, our calculations are relevant only at tempera- 
tures higher than the Kondo temperature Tk- 

The method can be applied to study other many-body 
effects that are expected to be important in nanoscopic 
structures due to strong electron-electron and electron- 
phonon coupling. A more general implementation of the 
presented method is under way. 
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